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INTRODUCTION 



in 
o 



A dirac equation with an interaction linear in coordinates was considered long ago 1] and recently rediscovered in the 
context of the relativistic many body theoriesQ. The equation is named Dirac oscillator, since in the non-relativistic 
limit it becomes a harmonic oscillator with a very strong spin-orbit coupling term. Dirac oscillator has attracted much 
attention and the concept gave rise to a large number of papers concerned with its different aspects |3t 14. l5l 1(1 1 71 1 81 1 9. 1 1 Cj . 

During the last years two-dimensional electron systems have become an active research subject due to the rapid 
growth in nanofabrication technology that has made possible to the production of low dimensional structures like 
quantum wells, quantum wires, quantum dots etc.[lllll3.ll3j. In non-relativistic case, the two dimensional parabolic 
potential has often been used to describe the spectrum of the electron in confined two-dimensional systems. For the 
relativistic case, the spectrum and properties of the such systems can be determined by using two dimensional Dirac 
oscillator[13|lJ,|l5|. Despite their simplicity, both relativistic and non-relativistic oscillator potentials appear to be 
a good approximation to complicated low dimensional nanostructures. 

In order to present a more realistic model we construct deformed form of the Dirac oscillator including a term 
qr 3 . The Dirac oscillator with a deformed term becomes radial sextic oscillator potential in the non-relativistic limit. 
The non-relativistic sextic oscillator potential is quasi-exactly-solvable(QES) for which it is possible to determine 
algebrically a part of spectrum but not whole spectrum [lil UtI Hli ] . It will be shown that the solution of the relativistic 
^ | ■ sextic oscilator can also be treated in the context of the QES problem. 

We also investigate the effect of the magnetic field on the Dirac sextic oscillator. It will be shown that for a specific 
values of the magnetic field the Dirac sextic oscillator problem is exactly solvable. 

The paper is organized as follows. In section 2 we discuss the construction of the Dirac sextic oscillator in polar 
coordinate, in the presence of the magnetic field. In section 3, we solve the Dirac sextic oscillator in the context QES 
problem. We show that the problem possesses hidden s^— algebra. In section 4 we briefly conclude our method and 
results. 

On 

construction of the dirac sextic oscillator 

in 

The (2 + 1)— dimensional Dirac equation for free particle of mass M in terms of two-component spinors ip, can be 
written as 

-i— > 
cr. 



i> (i) 



Since we are using only two component spinors, the matrices /3 and (3ji are conveniently defined in terms of the Pauli 
spin matrices which satisfy the relation 0%(Tj = Sij + SijkO~k, given by 



071= (7i; /372=<r 2 ; = a 3 . (2) 

In (2 + 1)— dimensions, the momentum operator pi is two component differential operator, p = —ih(d x ,d y ) for free 
particle. In the presence of the magnetic field it is replaced by p — > p— eA, where A is the vector potential, and 
the 2D Dirac oscillator can be constructed by changing the momentum p — > p— imuj(3r. We are now looking for 
some expressions on the right hand side of Q that can be interrepted as radial sextic oscillator Hamiltonian in the 
non-relativistic limit. For this purpose we introduce the following momentum operator 

p — > p— icr 3 (mu; — qr 2 )r (3) 

where q is a constant and the Dirac equation takes the form 

[E - a mc 2 ] ip = c [cr + (p x - ip y - i(Mu) - qr 2 )(x - iy))] ijj + 

c [er_ (p x + ip y - i(Muj - qr 2 )(x + iy))] tfj (4) 

In polar coordinate, x = rcosifi, y — rcos<j), the 2D Dirac equation can be written as 

d 2 Id Id 2 2t (Afw - qr 2 ) d 



e 2 ^ = - c 2 h 2 



4>i 



dr 2 r dr r 2 dcf> 2 fi d<f> 

+c 2 (q 2 r 6 - 2Mcjqr 4 + (M 2 uj 2 - ATiq) r 2 + 2hMu) ifa (5) 
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where if>i = ipx(r, (f>) is the upper component of the spinor ip and e 2 = E 2 — M 2 c 4 . The substitution of 

p—iTntfi 

^(r,0) = -^-/(r) (6) 

leads to the following equation 



2t2 9 2 /(r) 



where V(r) is given by 



e z f( r ) = -c'r-^ + [V(r) + 2ftMw(l - m)] /(r) (7) 



ft 2 f m 2 _ I) 

V(r) = — — — + q 2 r 6 - 2Mujqr A + (M 2 uj 2 - 2hq{2 - m)) r 2 (8) 



It is not difficult to see that in the non-relativistic limit J7J corresponds to the Schrodinger equation with radial 
sextic oscillator potential with a spin-orbit coupling term. In the presence of the symmetric gauge vector potential 
A(x,y) = B/2(—y,x,0), the potential (JSJ takes the form 

n 2 (m 2 - i) 
V(r) = — K — — ^- + fieB{m-l) + 

q 2 r 6 - {2Mu - eB) qr 4 + ( (mu) - - 2%q(2 - m)J r 2 (9) 

From now on we restrict ourselves to the solution of Q and (jHJ. 

METHOD OF SOLUTION 

In this section we show that the Dirac sextic oscillator (JJJ is one of the recently discovered quasi-exactly solvable 
operator 0, 0] . It is well known that the underlying idea behind the quasi exact solvability is the existence of a 
hidden algebraic structure. Let us introduce the following realization of the sl2— algebra: 

dp dp dp I 

The generators satisfy the commutation relations of the sli— algebra for any value of the parameter j. If j is a positive 
integer the algebra (|l(Jfl possesses j + 1— dimensional irreducible representation: 

P j+1 = (I, p, p?) (11) 

The linear and bilinear combinations of the operators given in (|10fl are quasi exactly solvable, when the space is defined 
in (|llfl • In order to show the Dirac sextic oscillator has a s^-symmetry, let us consider the following combinations of 
the operators ljTU|) : 

T = -J J- + J ~^r L J- + 16c 4 h 3 qJ + + Ac 2 TiMluJ . (12) 

Then, the eigenvalue problem can be written as 

TP k {p) = (e 2 + 2Mc 2 huj) P k (p) (13) 

where Pk (p) is the k th degree polynomial in p. Let us turn our attention to the Dirac sextic oscillator J2J • Introducing 
a new function 

/(r) = r m ~h-^ r2 -^F(r) (14) 
and then changing the variable r = 2chy/~p, we obtain the following expression 

e 2 F(p) = 

-p^M + ( m -l + Ac 2 hcoMp - 16 9 cW) 

op 1 op 

-c 2 (4Mhuj(m - 1) - 16qc 2 h 3 (m - 2)p) F(p). (15) 
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We can show that the eigenvalue equation J7J and 113(1 are identical, when the following holds: 

m=i + 2, F(p) = P k (p). (16) 
When the generators act on the polynomial l|ll|) , we can obtain the following recurrence relation 

16gc 4 ?L 3 (fc - j)P k +i(e) + (e 2 + 4Mc 2 huj(j -k+ l))i\(e) - 

fc(i-A; + 2)P fe _ 1 (A) =0 (17) 

with the initial condition Po(l) = 1- If £ i is a root of the polynomial Pfc+i(e), the wavefunction is truncated at 
k = j and belongs to the spectrum of the Hamiltonian T. This property implies that the wavefunction is itself the 
generating function of the energy polynomials. The roots of the recurrence relation l|17l) can be computed and the 
first few of them are given by 

Pi(e) = e 2 + AMc 2 hu 

P 2 (e) = ( £ 2 + AMc 2 huj) (e 2 + 8Mc 2 Huj) - 32qc 4 h 3 

P 3 (e) = e 6 + 8c 2 h(3e 4 Muj + 48c 4 h 2 Muj{M 2 uj 2 - 3qh) + (18) 

2c 2 he 2 (11 M 2 lu 2 - 10qh)) 
Pile) = e 8 + AQe & Mc 2 hoj + 8Qe 4 c 4 Ti 2 {7M 2 uj 2 -6Tiq) + 

128e 2 c 6 h 3 Muj(25M 2 Lj 2 - 69%) + 6144c 8 ?i 4 (3?iV - QM 2 uo 2 fiq + M 4 uo 4 ). 

The function Pj(p) forms a basis for sl2— algebra and it can be written in the form 

p 3 (p)=iz ckP ^ £ )p k - ( i9 ) 

fc=0 

In the presence of the magnetic field the Hamiltonian can be solved by the same procedure given above. When the 
magnetic field B = 2Mu/e, the Hamiltonian (JJJ) takes the form 

■VCD--**'™- 



c 2 



dr 2 

+ q 2 r 6 - 2hq(2 - m)r 2 + AHMuj(1 - m) 



f(r) (20) 



We define the wavefunction 



f(r) = r*~ m e- ^ F{r) (21) 
and changing the variable r = 2ch^fp, we obtain the following differential equation 

e 2 F{ P ) = -P^ + (i + 1 - 16qc 4 h 3 p 2 ) + lQ jqc ^pF{p). 

When F(p) = Pk(p), where Pk{p) is k th degree polynomial in p, we obtain the following reccurence relation 

16gc 4 ?L 3 P fc+1 (£) - e 2 P k (e) + k(j + 2 - fc)P fe _ 1 ( £ ) = 

The polynomial Pk(e) vanishes for k = j + 1 and the roots of Pj(e) belongs to the spectrum of the l(2()|l . We list the 
first ten of them below 





= e 2 




P2{e) 


= e 4 - 


- 2rf 




= £ 6 - 


- 10eV 


P 4 (e) 


= £ 8 - 


- 30eV + 72t7 4 




= e w 


- 70eV + 712eV 


Pe(e) 


= e 12 


- 140e 8 77 2 + 382(teV - 1080077 6 


Pr(e) 


= e 14 


- 252e 1( V + 14796eV - 164592eV 


Ps(e) 


= e 16 


- 420c 1 V + 46380eV - 1307600eV + 4233600?7 8 


P S (e) 


= e 18 


- 660e 14 ?7 2 + 125004e 1 V - 7250320eV + 88504707e 2 ?7 



(22) 
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where rj 2 = 16gc 4 ?i 3 . Therefore we can obtain the eigenfunction of i|20[l in the closed form. We conclude that 
anhormonic interaction destroys the general symmetry of the Dirac equation, but the specified magnetic field can 
restore the symmetries of the Dirac equation. This feature implies that for the specific values of the magnetic field 
B = 2hu/e, analytical solutions of the roots of the polynomials are available. 



CONCLUSION 



Due to the interest of the lower dimensional field theory and condensed matter physics, we have constructed Dirac 
sextic oscillator in two-dimensional space in the presence of the magnetic field. We have given eigenstates of the 
corresponding equation interms of the orthogonal polynomials ^IJ. It has been shown that the Dirac sextic oscillator 
possesses hidden sl2— symmetry. 

We note that similar constructions of the other quasi-exactly solvable Dirac equations also seem possible by consid- 
ering the momentum operator of the form: p — > p— ia^(muj + v(r))r. The method we have introduced can be applied to 
construct Dirac equation including the other non-relativistic potentials such as Poschle- Teller, Eckart, Scarf, Hulten, 
etc. We hope that the Dirac sextic oscillator may be used as a model in the related fields of the physics. 
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